sigma_000 = sqrt(kappa1)*(alpha(:,1,1,1) + beta(:,1,1,1));
sigma_011 = sqrt(kappa1)*(alpha(:,1,2,2) + beta(:,1,2,2));
sigma_001 = sqrt(kappa1)*(alpha(:,1,1,2) + beta(:,1,1,2));
sigma_111 = sqrt(kappa1)*(alpha(:,2,2,2) + beta(:,2,2,2));
hold on;
plot(sigma_000,'r','LineWidth',2); plot(sigma_011,'r--','LineWidth',2); plot(sigma_111,'b','LineWidth',2); plot(sigma_001,'b--','LineWidth',2); 
